On the connection between shape and stellar population in 

early-type galaxies 



M. D'Onofrio 1 and T. Valentinuzzi 1 

Astronomy Department, Vicolo Osservatorio 3, 1-35122 Padova, Italy 
mauro . donof rioOunipd. it 

G. Fasano 2 and A. Moretti 2 and D. Bettoni 2 and B. Poggianti 2 and B. Vulcani 1 ' 2 and J. 

Varela 2 

IN AF/ Astronomical Observatory of Padova, Vicolo Osservatorio 5, 1-35122 Padova, Italy 

J. Fritz 3 

Sterrenkundig Observatorium, University of Gent Krijgslaan 281 S9, B-9000 Gent, Belgium 

A. Cava 4 

Instituto de Astrofisica de Canarias, via Lactea E38205, La Laguna (Tenerife) , Spain 

P. Kjaergaard 5 

The Niels Bohr Institute for Astronomy Physics and Geophysics, Juliane Maries Vej 30, 2100, 

Copenhagen, Denmark 

M. Moles 6 

Centro de Estudios de Fisica del Cosmos de Aragon (CEFCA), G/ General Pizarro, 1, 44001 Teruel, Spain 

W.J. Couch 7 

Center for Astrophysics and Supercomputing, Swinburne University of Technology, PO Box 218, Hawthorn 

Victoria 3122, Australia 
and 

A. Dressier 8 

Observatories of the Carnegie Institution, 813 Santa Barbara Street, Pasadena, California 91101, USA 



ABSTRACT 

We report on the discovery of a relation between the stellar mass M* of early- type galaxies 
(hereafter ETGs), their shape, as parametrized by the Sersic index n, and their stellar mass- 
to-light ratio M* /L. In a 3D log space defined by these variables the ETGs populate a plane 
surface with small scatter. This relation tells us that galaxy shape and stellar population are 
not independent physical variables, a result that must be accounted for by theories of galaxy 
formation and evolution. 

Subject headings: galaxies: general, galaxies: fundamental parameters, galaxies: structure 
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1. Introduction 

It has long been known that the shape of galax- 
ies is c onnected with t heir stellar population prop- 
erties. iMorgan (119581) based his galaxy classifica- 
tion scheme of spiral galaxies on the relation be- 
tween the central light concentration and spectr al 
characteristics (see also Morgan fe Mavalllll957l ). 

ETGs have been considered a homogeneous 
class within the Hubble sequence, manifesting 
a restricted range of shapes and a generic old 
stellar population. Their internal structure is 
usually described as differing by only a simple 
scale factor from galaxy to galaxy (structural ho- 
mology). Observational support for this inter- 
pretation was based on the fact that all lumi- 
nosity profiles appeared to follow the de Vau- 
coule urs' law (de Vaucouleurd 1948 ) I(r) ~ r 1 / 4 . 
First, lMichardl(ll98l " and lSchombertl (| 19861 ) noted 
that residuals of the individual r 1 / 4 fits corre- 
lated with galaxy luminosity. Then, the first 
breakthrough in the previously outlined scenario 
ca me more tha n 20 ye ar s ago through the papers 
by Capacciolil (1987); de Carvalho fc da Costa 



(119881 ); ICapacciolil ([19891 ); iBurkertl ([1993I ). The 
final proof of structural non-homology came when 
better fits of the luminosity profiles were obtaine d 
using a Sersic' law I(r) ~ r 1 /™ (jSersid Il968l) . 
where n is a free index that defines the shape of the 
light profile and correlates with the to tal luminos- 
ity of the galaxy: the n — L r elation flCaon et al 



Young fc Currid Il994t iPrugniel fc Simien 



1993; 
1997). 

It is also well known that ETGs share a 
3D planar distribution with a small scatter, 
the so-called Fundamenta l Plane (hereaf t er FP , 



Diorgovski fc Davisl 119871; iDressler et all 119871) . 
defined by the effective radius R ei the effective sur- 
face brightness (^) e , and the central velocity dis- 
persion a. The observed tilt of the FP relative to 
the virial theorem expectation, and the small scat- 
ter around the plane, are still matters of debate. 
Many researchers believe that the tilt is likely 
a consequence of stellar population effects which 
manifest themselves through the observational ev- 
idence that low-mass ETGs tend to be younger 
and to show lower stellar mass-to-light ratios with 
respect to the massive ones (the M* /L — L re- 



lation; see e.g. Faber et all 119871 iGerhard et al 



20011 iTreu et al.H2005l ). Actually, the observa- 



tions reveal that the tilt is still substantial in the 
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see e.g. 



Pahre et al 



19981 ). where the effect of the stellar mass-to-light 
ratio variation should be mild. Another possible 
explanation of the tilt is linked to the amount and 
distribution of dark matt e r (DM ). According to 
ICiotti. Lanzoni fc Renzinil (|l996T ) the DM expla- 
nation requires a peculiar fine-tuning to reproduce 
the observed p roper ties of the FP. More recently 



Tortora et al 



(|2009l) . estimating the total M/L 
from simple Jeans dynamical models, found that 
an important role should be played by DM, since 
its fraction within R e seems to be roughly constant 
for galaxies fainter than Mb ~ —20.5, while it in- 



creases for brighter galaxies. iPadmanabhan et al.. 
(|2004l) and iHvde fc Bernardil (|2009l) also found 



evidence that the M tot /M* ratio increases with 
mass, a fact that may imply a dark-to-bright mat- 
ter ratio increasing along the FP. 

Another popular explanation of the FP tilt, 
alternative (or complementary) to the previous 
ones, rests upon the already mentioned evidence of 
structural non-homo logy of ETGs, parametrized 
by the Sersic index n dBertin. Ciotti fc Del Principd 



20021 iTruiillo. Burkert. fc Bellll2004l) ~ 

Here we report one of the re sults of a work in 
progress (|D'Onofrio et al.ll2010l) . indicating that a 
connection between galaxy shape and ste llar pop- 
ulation that was found by[M 

organ (1958) for spi- 
ral galaxies, is also present in ETGs. We suggest 
that such a relationship is likely connected to the 
properties of the FP, in particular with its scatter, 
and is therefore one of the fundamental physical 
relations that underpins its origin. 

2. Our galaxy sample selection 

The present work is based on the data of 
the Wide-field Nearby Galaxy-clusters Survey 



(WINGS, see e.g. 1 


? asano et al. 2006: Varela et al. 


2009: Cava et al. 


2009: Valcntinuzzi et al.l 2009: 


D'Onofrio et al. 20081). Our sample comes from 



the cross-match of the V and K band surface 
photometry. Total luminosities, effective radii and 
Sersic indices were measured with the automatic 



softw are GASPHOT (jPignatelli. Fasano fc Cassata 
20061 ). The quality of the GASPHOT photometry 
(magnitudes, radii, surface brightness and flatten- 
i ng) is discussed in detail in several WINGS papers 
(|Fasano et al.ll2006l : fVarela et alJliool ICava et aT 
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2001 IValentinuzzi et all l2009t 
20081; iFritz et al.ll2007ll2010h . 



D'Qnofrio et al 



Ga laxy masses have been derived by lFritz et al 



( 201dh fitting the observed (optical ) spectra of the 
WINGS survey (|Cava et all 120091) with theoret- 



ical spectra, obtained summing up Simple Stel- 
lar Populations (SSP) spectra of different ages, 
each one proper ly extinguished a ccording to the 
SSP's age itself (jFritz et al.ll2007l) . Details about 
the models and the fitting procedure (evolution- 
ary tracks, mass range and evolutiona ry stages of 



stars, IMF, SFH, etc.) can be found in lFritz et al 



(|2007Ll2010h . Here we just mention that the equiv- 
alent widths of the main lines and the flux emit- 
ted in significant ranges of the spectral contin- 
uum have been used to find the best fit of the 
observed spectra, at varying the mass fractions of 
the SSPs, as well as the "selective extinction" and 
the metallicity. The masses obtained in this way 
have been found in fair agreement with those ob- 
tained from our GA SPHOT photometry, fo llowing 
the prescriptions of Bell fc de Jond ( 2002 ). 

All galaxies were classified as early-types (pure 
ellipticals or SOs) using MORPHOT, an automatic 
morphological classification pa ckage purposely de- 
vised for the WINGS pr oject (jFasano fe Vanzeha 
20071: iFasano et al.ll2010t) . 

Galaxy magnitudes, effective radii and Sersic 
indices have been fo und in good agreement both 
with the SDSS data (jValentinuzzi et alj|2010f ) and 
with the values w e obtained (for co mparison) us- 
i ng the GALFIT dPeng et alj |2002l) and GIM2D 
( Simard et al. 20021 ) packages. 

From an initial sample of 2613 (848) galaxies 
with V (K) band data, we extracted a second 
" high quality sample" consisting of 240 ETGs with 
good photometry in both bands, well defined mor- 
phological types, optimum Sersic light profile fits 
and robust stellar mass measurements. The se- 
lected galaxies are all spectroscopically confirmed 
cluster members. 

3. The log(M*) - log(n) - log(M */L) relation 

The n — L and M*/L — L scaling relations 
mentioned in the introduction have been extcn- 



PniR 


niel & Simienlll997l Ciottil 1991 : Busarello et al.l 


1997 


Truiillo. Burkert. & Belli 2004| Grab 


am & Collcss 


1997 


Bertin. Ciotti & Del Principel2002: 


Faber et al. 



19871 iGerhard et al.ll200lHTreu et al.ll2005l) . They 
provide evidences that faint ETGs are character- 
ized by small values of n and M* /L with respect 
to bright objects. The origin of such behaviour is 
still a matter of debate. 

In this letter we highlight the connection be- 
tween n (i.e. the shape of the light profiles) and 
M* /L (i.e. the stellar populations) which is hid- 
den in the M* - n- M*/L relation. The link be- 
tween these two fundamental variables that char- 
acterize ETGs was indeed poorly analyzed in the 
literature, since the two variables do not correlate 
directly. 

The upper panels of Fig. [T] show the log(n) — 
log(M*) and \og{M* / L) - log(M*) relations for 
the "high quality sample". The lines mark the 
bilinear least-square fits to the data obtained us - 
ing the SLOPE program (jFeigelson fc Babulll992t ). 
The lower panels of the same figure show the 
orthogonal residuals of the two relations versus 
log(M*/L) and log(n), respectively. Note that the 
same trend of the residuals is observed when the 
classical least-square fit is adopted and the residu- 
als are evaluated along the Y axis. It is clear from 
the plots that the residuals of the log(M*/L) — 
log(M*) and log(ra) - log(M*) relations depend, 
respectively, upon log(n) and log(M*/L). These 
two variables are independent of each other, but 
are connected through their residuals. The same 
mass-to-light ratio is achieved by galaxies with 
a low/high mass when their Sersic index (n) is 
high/low, while the same value of n is achieved 
by galaxies of low/high masses when M* / L is 
high/low. Despite the well known difficulties in 
measuring the Sersic index (variance An ~ 1), the 
observed trends reveal a clear interplay between 
galaxy shape, mass and mass-to-light ratio. This 
behavior of the residuals implies that ETGs are 
distributed on a plane in a 3D space whose axes 
are log(M*), log(rc) and \og{M*/L) (see Fig. EJ. 
This relationship is observed in both V and K 
bands and for both the " high quality" and the to- 
tal samples. In particular, for the "high quality" 
sample the equations of the planes which minimize 
the sum of the distance from the fit are: 



log(M*) = 1.12(±0.22)log(M7£) v 
+2.19(±0.23) log(rv) 
+9.14(±0.11) (1) 
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Fig. 1. — (Upper panels) The log(n) — log(M*) 
(left) and \og(M* / L) - log(M*) (right) relations. 
The gray scale of the dots mark low(white) to 
high(black) values of log(M*/L) and log(n) re- 
spectively. The lines show the best fits. (Lower 
panels) The orthogonal residuals of the above re- 
lations plotted versus log(M*/L) and log(n) re- 
spectively. 





Fig. 2.— The plane formed by ETGs in the 3D 
space defined by log(M*), log(n) and log(A/*/i) 
seen along two different viewing angles. 



log(AT) = 1.77(±0.55)log(M7i) K 
+2.11(±0.26) \og(n K ) 
+9.70(±0.14), (2) 

for the V (rms = 0.10) and K (rms = 0.11) 
bands, respectively. The Principal Component 
Analysis (PCA) confirms that we have a planar 
distribution. The first two eigenvectors for the V 
band are: 

£ = 0.67 log(M*) + 0.59 log(rc) 

+0.45 \og(M*/L) (3) 

X] = 0.08 log(M*) +0.54 log(ra) 

-0.831og(M*/L). (4) 

They are built with the three variables, and to- 
gether support up to ~ 90% of the total variance 
(respectively ~ 60% and ~ 30%). The same is 
true for the K band with small differences in the 
PCA coefficients. 

As far as we know, the n — M* relation for 
ETGs is explicitly presented here for the first time, 
at least for ETGs. Actually, iDuttonl (|2009t ) pub- 
lished the n — M* relation for disk galaxies in the 
redshift range 0.005 < z < 0.05 using the data 
of the New Yo rk University Value Added Catalog 
fNYU-VAGC. lBlanton et~aDl2005h . Their slope is 
a bit lower than our, likely because of the limited 
range of values of n for disk objects (from 0.2 to 0.6 
in log scale). However, remarkably enough, their 
Fig. 5 clearly shows that the residuals with respect 
to a straight line fit depend on colors, i.e. on the 
M*/L ratio. This suggest that the M*-n-M*/L 
relation might be valid over the whole Hubble se- 
quence. 

We can comp are our n — L and n — R e rela- 
tions with that of Caon et al. ( 19931) , based on the 
Virgo and Fornax ETGs. For the V band we have: 
log(n) = 0.60(±0.05)log(L) - 5.74(±0.3) with a 
correlation coefficient C.C. = 0.51 and log(i? e ) = 
0.73(±0.06)log(n) + 0.02(±0.03) with a C.C. = 
0.54. Both relations are in agreement (within 3a) 
with t he previous determinations by ICaon et al 
(Il993h . Actually, in the case of the n — L relation 



our slope is larger than in lCaon et al.l (|1993l 0.6 vs. 
0.45). However, we must consider that the slope 
value depends on the adopted fitting strategy. In 
fact, using a simple least-square fit, we get the re- 
lation log(n) = 0.30(±0.05)log(L) - 2.42(±0.3), 
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quite different from that o btained when a bili near 
LSQ is used. The data of ICaon et al.l (|l993l ) are 



based on 40 ETGs against the 240 of the WINGS 
database. Since the errors are approximately of 
the same order for both variables, we believe that 
the bilinear LSQ fit is more appropriate for deter- 
mining the correct slope of the relation. Therefore 
we estimate that the true slope of the n — L rela- 
tion is in the interval 0.4 -4 0.6. 

In the V band our M*/L — M* relation has 
a slope of ~ 0.6 (again the slope depends on the 
fitting strategy adopted, varying from ~ 0.45 (for 
the classic LSQ fit) to ~ 0.6 (for the bilinear fit)), 
while the M* /L — L relation has a slope of ~ 0.03 
(with a C.C.~ 0. 03), in agreement wit h the recent 
determination bv lTortora et al. I (|2009h . Note that, 
given the slope of our M*/L - M* relation (0.6), 
some trivial algebra would imply that M* / L ~ L 2 , 
largely at variance with the formal value we found 
(2.0 vs. 0.03). However, it is easy to recognize 
the origin of this apparent discrepancy in the lack 
of correlation between M*/L and L (C.C.~ 0.03), 
which makes practically undefined the slope of the 
best fitting straight line. Such lack of correlation 
is also the reason why there is a 3D planar dis- 
tribution of ETGs in the M * - n - M*/L log 
space, while it is not present in the correspond- 
ing L — n — M*/L space. 

Since we were not able to find in the litera- 
ture any estimate of the slope of the M* /L — M* 
relation, we tried to get an 'indirect' check of 
our value combining the slopes of the a — M* 
and M* /L — a relations. Our data give a ~ 
o.26-=-o.38 ^ depending on the adopted fitting 
tool (thes e values are within th e range of slopes 



given by iGraves fc Faberj l2010t) . and M*/L ~ 
cr 4 ^ 1 - 3 . The combination of these equations pro- 
vides M*/L ~ j^*o.i-r-o.6 ) m a g reem ent with our 
'direct' finding. In any case, it is worth stress- 
ing again that the existence of the 3D M* — n — 
M* /L relation does not depend on the final slopes 
adopted for the M* - n and M*/L - M* relations. 

4. Consequences for theories of galaxy for- 
mation and evolution 



Be fore the discove r y of th e FP ( Diorgovski fc Davis! 



1987b 

and (u) P . 



Dressier et al.l 11987), the L — a, a — R, 



of ETGs. When lTerlevich et al.l (|l98lh found that 
the residuals of the Faber-Jackson relation were 
connected with the effective surface brightness 
(fi) e , the existence of a 3D relationship between 
(j—R e — {n) e (i.e. the FP) was soon revealed. Sim- 
ilarly, we report here on the existence of a plane 
preferentially shared by ETGs in the 3D space 
defined by the quantities M*, n and M* /L, which 
are already known to be linked by mutual indepen- 
dent correlations. We define this plane and claim 
that its very existence implies that stellar popu- 
lations and shapes of ETGs are not independent 
quantities. Why this occurs is presently unknown, 
but this mutu al dependence recalls the correla- 
tion found by iMorgan ( 1958 ) for spiral galaxies 
and could be an important ingredient for theo- 
ries of galaxy formation and evolution. A 'nat- 
ural' guess one can make is that such a relation 
is connected to the FP of ETGs, i.e. with the 
process of virialization that produced, after the 
gas collapse phase, the final stellar structure in 
dynamical equilibrium. How can this finding be 
reconciled with the hierarchical merging scenario 
of galaxy formation where ETGs are the final out- 
put of a series of merging events? How is this 
relationship connected with the well known cor- 
relations of galaxy mass and mass-to-light ratio 
with age and metallicity? How it is connected 
with the FP? We are presently investigating the 
implications of our finding. Preliminary analy- 
ses based on FP simulations in which the Sersic 
indices of the mock galaxies randomly share the 
observed distribution, or (alternatively) are forced 
to obey (with the proper scatter) the relation (1), 
suggest that it might act as a sort of 'fine-tuning' 
effect that keeps small the scatter around the FP 
(jD'Onofrio et al.ll2010h . 



Benedetta Vulcani and Bianca Maria Poggianti 
acknowledge the financial support from the ASI 
contract 1/016/07/0. We also thank the anony- 
mous referee for his help in improving this work 
a nd for br i nging our attention toward the paper 
of lDuttonl (|2009l ). 



R P relations (jFaber fc Jackson 1976; 



Kormendvlll977h were well known scaling relations 
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